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Abstract 

We investigate flavor violation mediated by Higgs bosons and supersymmetric particles in a 
predictive class of models based on the non-Abelian flavor symmetry Qq. These models, which 
aim to reduce the number of parameters of the fermion sector and to solve the flavor changing 
problems of generic SUSY setup, assume three families of Higgs bosons and spontaneous/soft 
violation of CP symmetry. Tree-level contributions to meson-antimeson mixings mediated 
by Higgs bosons are shown to be within experimental limits for Higgs masses in the (1-5) TeV 
range. Calculable flavor violation induced by SUSY loops are analyzed for meson mixing and 
lepton decays and found to be consistent with data. Significant new SUSY contributions arise 
in B s — B s mixing, but non-standard CP violation is suppressed. A simple solution to the 
SUSY CP problem is found, which requires light Higgsinos. 
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1 Introduction 



The gauge interactions of the standard model (SM) fermions are invariant under separate U(3)l x 
U (3)# transformations. This global symmetry is broken explicitly by the fermion Yukawa couplings. 
For the light fermions violation of this symmetry is small, being proportional to their masses. This 
feature has played a crucial role in the success of the SM in the flavor sector. In extensions of 
the SM this property is generally lost, often leading to excessive flavor changing neutral current 
(FCNC) processes. 

A case in point is the supersymmetric standard model which is the subject of this paper. While 
the gauge interactions of the SUSY SM respect the U(3)l x U(3)r global symmetry, there are 
new sources of violation of this symmetry, in the soft SUSY breaking sector. Indeed, generic soft 
SUSY breaking scenarios lead to excessive FCNC in processes such ad K° — K° mixing, B° — B° 
mixing, D° — D° mixing, and flavor changing leptonic decays such as fi — > &-/ [1]. This problem 
is most severe in the K° — K° system. SUSY box diagrams involving gluino and squarks modify 
the successful SM prediction for AM^ and e#, leading to the following constraints for the real and 
imaginary parts of the amplitude [2]: 

| (Re, Im)(^ i ) 12 (4 ii ) 12 | 1/2 < (9.6 • 1(T 4 , 1.3 • 1(T 4 ) C^^j • (1.1) 

Here (5ab)i? = (^abIw'M 2 * s a fl avor violating squark mass insertion parameter, for (A,B) = 
(L, R), with m being the average mass of the relevant squarks (d and s in this case). For this esti- 
mate the gluino mass was assumed to equal the average squark mass. Now, the natural magnitude 
of the mixing parameters {Sf^u and (5^) 12, in the absence of additional symmetries, should be 
of order the Cabibbo angle, ~ 0.2. Since the parameters (Sab)^ split the masses of the squarks, 
one sees from Eq. (II. ip that a high degree of squark mass degeneracy is needed for consistency. 
Analogous limits from B® — B® mixing are less severe, as given by [3] : 

I (Re, Im)(^) 13 (4*)i3| 1/2 < (2-1 • 10" 2 , 9.0 • lO" 3 ) (^^7) ■ (1-2) 

Note that the natural value of this mixing parameter, in the absence of other symmetries, is 
V u b ~ 3 x 10~ 3 . The constraints from Eq. (11.21) are well within limits. B s — B s mixing provides 
even weaker constraints. 

It can be argued that a natural explanation for solving this problem is to enhance the symmetry 
of the SUSY SM by assuming a non-Abelian symmetry G (a subgroup of the U{3)l x U{3)r\ 
that pairs the first two families into a doublet, with the third family transforming trivially [5]Jj 

3 Grouping all three families into an irreducible triplet representation of G is also possible. The large top quark 
mass however reduces the original U(3)l x U(3)r symmetry to U(2)l x U(2)r, so we find it is easier to work with 
(2 + 1) assignment. 
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Invariance under G will then lead to degeneracy of squarks, as needed for phenomenology. A variety 
of such models have been proposed in the literature [5], [6], [7], [8], [9]. In Ref. [5], SU(2) family 
symmetry and its variants were proposed to solve the SUSY FCNC problem. If the symmetry is 
global, one has to deal with the Goldstone bosons associated with its spontaneous breaking. Global 
symmetries are susceptible to violations from quantum gravity. Local gauge symmetries are more 
natural, but in the SUSY context there would be new FCNC processes arising from the family 
SU(2) .D-terms [10]. Exceptions to this generic problem are known to exist [6]. 

A more natural solution to the problem is perhaps to choose G to be a non-Abelian discrete 
symmetry group [7]. In this case there would be no Z)-term problem, since there are no gauge 
bosons associated with G. Spontaneous breaking of such symmetries will not lead to Goldstone 
bosons. If the symmetry breaking occurs before the inflationary era, such models should also be 
safe from potential cosmological domain wall problems. Such non-Abelian discrete symmetries 
have found application in understanding the various puzzles associated with the quark and lepton 
masses and mixing angles with or without supersymmetry [11], more recently for understanding 
the tri-bimaximal neutrino mixing pattern [12]. It would be desirable to find a symmetry that 
sheds light on the fermion mass and mixing puzzle, and at the same time solves the SUSY FCNC 
problem. 

The supersymmetric standard model has another problem. In the flavor conserving sector CP 
violation is generically too large. Neutron and electron electric dipole moments (EDM) receive 
new contributions from SUSY loops. Unless the new phases in the SUSY breaking sector are small 
or conspire to be small, experimental limits on the EDM of the neutron (d n ), electron (d e ), and 
atoms will be violated by two to three orders of magnitude (depending on the squark and slepton 
masses) [13], [14]. The imaginary parts of the left-right squark mixing parameters must satisfy the 
constraints (from the experimental constraints d n < 6.3 x 1CT 26 e-cm, d e < 4.3 x 1CT 27 e-cm) [15] 

M( < j 8 ) I j<i.9xio-( 5j ^), MfJU.Jsi.rxiir'^), (i.3) 

assuming that the gluino/Bino has the same mass as the squark/slepton. Now, since these mixing 
parameters are expected to be suppressed by fermion helicity factors (but enhanced by the MSSM 
parameter tan/5) the natural values for these mixing parameters are of order (1 x 10~ 4 , 3 x 1(U 6 ) 
respectively, (for tan (3 = 10 and assuming order one phases). This implies that the CP violating 
phases arising from the soft SUSY breaking sector must satisfy 9 d < 1/53, 9 e < 1/63 (for gluino 
(Bino) mass of 500 GeV (100 GeV). Why this is so, while the Kobayashi-Maskawa phase takes 
order one value, is the SUSY CP puzzle. It would be desirable to resolve this puzzle based on a 
symmetry principle in the same context where the SUSY FCNC problem is solved. 

The purpose of this paper is to study a recently proposed SUSY model based on the non-Abelian 
symmetry group Qq [8] which addresses these issues. Q 6 is a finite subgroup of SU (2) with twelve 
elements. Apart from providing a solution to the SUSY flavor problem, this class of models can 
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also constrain the quark masses and mixings. It was shown in Ref. [8] that with the assumption 
of spontaneous (or soft) CP violation, there is a non-trivial relation between quark masses and 
mixings in this model. This sum rule was found to be consistent with experimental data. 

A crucial aspect of the Qq model relevant for the quark mixing sum rule is that CP violation 
occurs either spontaneously or softly. This can help ameliorate the SUSY CP problem mentioned 
above. CP invariance requires that the gaugino masses, the /i terms and the trilinear A terms 
be all real. In the Qq model of Ref. [8] it was found that there is a phase alignment mechanism 
that makes the phases of the sfermion mixing terms arising from the v4-terms to align with the 
phases of the fermion masses. So SUSY CP violation is suppressed to a large extent. However, 
spontaneously induced complex VEVs do lead to non-zero contributions to EDM. Here we analyze 
these contributions. Since these complex VEVs are accompanied by the Higgsino fi terms, a simple 
solution to the problem is found by making the Higgsinos to be lighter than the squarks. Adequate 
suppression of EDM is obtained for /i ~ 100 GeV, while squark masses are of order 500 GeV. This 
suggestion obviously has testable implications for physics that will be probed at the LHC. 

The fermion mass matrices that allow for a non-trivial prediction and the phase alignment is 
a generalization of well studied models [16]. The mass matrices for up and down quarks and the 
charged leptons take the following form: 

/ C \ 
M= ±C B . (1.4) 
\ B' A j 

The main feature of such mass matrices is that the phases can be factorized, i.e., M = P ■ M° ■ Q, 
with M° being real and P, Q being diagonal phase matrices. This feature, when combined with the 
Qq symmetry, has an the interesting consequence that CP violation induced by SUSY loops are 
suppressed. This will be discussed in more detail in Sec. 4. 

The form of Eq. (II. 4p can be obtained in renormalizable theories based on Qq symmetry. This 
requires the introduction of three families of Higgs doublets, which fall into 2 + 1 representations 
of the Qq group, very much like the quarks and leptons. With multiple Higgs fields coupling 
to fermions, invariably there will be tree-level FCNC mediated by the Higgs bosons. The flavor 
changing Higgs couplings are not arbitrary, but can be computed in terms of the fermion masses 
and mixings. We will show that these FCNC processes are within acceptable range, provided that 
the Higgs boson masses lie in the (1 — 5) TeV range (except of course for the standard model-like 
Higgs boson, which has a mass in the (100 — 130) GeV range). While Higgsinos are naturally light 
in this scenario, in the bosonic sector only the lightest SM-like Higgs will be accessible to LHC 
experiments. 

One of our major results is that non-standard CP violation is highly suppressed in this class 
of models. The phase factorizability of the fermion mass matrices implies that much of the SUSY 
induced CP violation is small. The structure of the Yukawa couplings in the model implies that the 
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amplitudes for tree-level FCNC induced by neutral Higgs bosons are nearly real (see discussions 
in Sec. 5). While there can be significant new contributions to meson-antimeson mixings, there is 
very little CP violation beyond the standard model. 

Our analysis is similar in spirit to that of Ref. [9]. Our approach is slightly different, with some 
differences in analytical results, fits, spectrum, and conclusions. In particular, we have presented 
complete analytical results for the Higgs boson spectrum, and we have a new proposal to solve the 
SUSY EDM problem, which requires light Higgsinos. We have also derived generalized constraints 
on SUSY FCNC parameters for the B^s — Bd,s system appropriate for a (2+1) mass spectrum. 

The plan of the paper is as follows. In Sec. 2 we describe the SUSY Qq model, lay out the 
parameter choice, and summarize the prediction for the quark sector. In Sec. 3 we analyze the 
Higgs potential involving the three pairs of Higgs doublets. We provide analytic expressions for the 
mass spectrum of Higgs bosons as well as numerical fits. Consistency of symmetry breaking and 
spontaneous CP violation will be established here. In Sec. 4 we address tree-level FCNC processes 
mediated by the heavy Higgs bosons. Sec. 5 is devoted to analysis of the SUSY flavor violation 
and EDM within the model. In Sec. 6 we conclude. 

2 Supersymmetric Qq Model 

Qq is the binary dihedral group, a subgroup of SU(2), of order 12. It has the presentation 

{A, B; A 6 = E,B 2 = A 3 , B~ X AB = A' 1 } . (2.5) 
The 12 elements of Q G can be represented as 

{E, A, A 2 , A 5 , B, BA, BA 2 , BA 5 }. (2.6) 
In the two dimensional representation the generators are given in a certain basis by 

/ cos ? sin ? \ „ / % \ , . 

A = -\ I ) B = n ■ ■ ( 2 - 7 ) 

\ — sin | cos | / V — % J 

The irreducible representation of Q$ fall into 2, 2', 1, l', 1 ", 1 ", where the 2 is complex-valued but 
pseudoreal, while the 2' is real valued. (Q 6 is the simplest group with two distinct doublet repre- 
sentations, which is very useful for model building.) The 1 and l' are real representations, while 
l" and l"' are complex conjugates to each other. The group multiplication rules are given as 

I' x 1' = 1, 1" x 1" = 1, 1'" x 1'" = 1', 1" x 1"' = 1, I' x 1'" = 1", 1x1" = I" 1 (2.8) 

2 x 1' = 2, 2 x 1" = 2', 2 x I'" = 2', 2' x 1' = 2', 2' x 1" = 2, 2' x 1"' = 2 (2.9) 
2 x 2 = 1 + 1' + 2', 2' x 2' = 1 + 1' + 2', 2x2' = l"+ Y" + 2 (2.10) 
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The Clebsch-Gordon coefficients for these multiplication can be found in Ref. [8]. 

The fermions of all sectors (up-quark, down-quark, charged leptons) are assigned to 2 + 1 
represtations of Qq. The model assumes three families of Higgs bosons, which are also assigned to 
2 + 1 under Qq. Their transformation properties are given by 

H=( I *}\=2', #3 = 1"'. (2.12) 



Here ip generically denotes the fermion fields, and H denotes the up-type and the down-type 
Higgs fields which are doublets of SU(2)l- Due to the constraints of supersymmetry, H u and H 3 
couple only to up quarks, while H d and H 3 couple to down-type quarks and leptons. The Yukawa 
couplings of the model in the down quark sector arise from the superpotential 

W = a d ^ c 3 H 3 + M T mP c 3 H - $ d ^ cT ir 2 H + 5 d ^ T mP c H 3 + h.c. (2.13) 

with similar results for up-type quarks and charged leptons. This leads to the mass matrix for the 
down quarks given by 

CO S d v d3 (3 d v d2 > 
S d v d3 (3 d v dl . (2.14) 
P'd V d2 P'd v dl UdVdZ ) 

Here v d i, v d2 , v d3 are the vacuum expectation values of Hf 2 3 fields, which break the Qq symmetry. 

Now, the potential of the Qq model admits an unbroken S 2 symmetry which interchanges 
H^ d *-> H 2 ,d . This unbroken symmetry allows us to choose a VEV pattern 

Vui = v u2 , v d i = v d2 , . (2.15) 

Consequently, a 45° rotation of the matrix in Eq. (12.141) in the 1-2 plane can be done both in the 
up and the down quark sectors without inducing CKM mixing. This will bring the mass matrices 
to the desired form of Eq. (11.41) . By using the unbroken 5*2 symmetry, we make a 45° rotation on 
the Higgs fields, H l<2 = (H 1 ± H 2 )/y/2, so that H x acquires a VEV, while (H 2 ) = 0. We shall drop 
the hat on these redefined fields, and simply denote the VEV of the redefined Hi as v\. 

We assume that CP is a good symmetry of the Lagrangian, and that it is broken spontaneously 
by the VEVs of scalar fields. If the full theory contains SM singlet Higgs fields, spontaneous 
CP violation in the singlet sector will show up as soft CP violation in the Higgs doublet sector. 
Explicit examples of this sort have been given in Ref. [8]. For now we simply assume that the 
Yukawa couplings in Eq. ( 12.131) are real, and the CKM CP violation has a spontaneous origin, via 
complex VEVs of the Higgs doublet fields. We denote the phase of these (redefined) VEVs as 

A6 U = arg(f M3 ) - &rg(v ul ), A9 d = axg(v d3 ) - axg(v d i). (2.16) 
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We make an overall 45° rotation on the Qq doublets, Q, D c and U c , and then a phase rotations 
on these fields: 

U -> P U U } U c -> P uC f/ c (2.17) 
and similarly for D and -D c fields, where 

/ 1 





P u , d = exp(z2A# Ujd ) 

V exp(iA0 u , d ) / 

/ exp(-i2A# nj d ) \ 

;c = 1 

\ exp(-iA0 u , d ) / 



(2.18) 



This will make the originally complex mass matrices of Eq. (11.41) real, which we parametrize as 



/ 



M u , d = m° t b 







lu, d/Vu, d \ 



-q u ,d/y u ,d o 
b' 



b u , d 

,2 



(2.19) 



V b ud yl d ) 

These real mass matrices can be diagonalized by the following orthogonal transformations: 





( 


m l, d 








\ 


Ol d M u , Ml, d O u , d = 
















\ 








m l 


J 




( 


m l, d 










Ou c , d c ^u, dMu c , d c O u c^ d c = 
















\ 








ml 


J 



(2.20) 



The CKM matrix Vckm is then given by 

VcKM = O^PqOd, 

where 



(2.21; 



(2.22) 



/ 1 \ 

e i2e " 
\ e w « ) 

with 6 q = A6 d - A9 U . 

Now it is clear how the Q§ setup reduces the number of parameters in the quark sector. The 
total number of parameters in the quark sector is nine (four real parameters each in M u and M d , 
plus a single phase 9 q ), which should fit ten observables. Spontaneous CP violation is crucial for 
this reduction of parameters. With explicit CP violation, there would have been one more phase 
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parameter. The single prediction of this model was numerically studied in Ref. [8], and shown to 
be fully consistent with data. Here we present a numerical fit to all the quark sector observables, 
which deviates somewhat from the fit given in Ref. [8]. The difference arises since here we have 
attempted to be consistent with the recent lattice determination of light quark masses. An excellent 
fit to the quark masses and mixings, including CKM CP violation, is obtained with the following 
choice of parameters at a momentum scale of /x = 1 TeV. 



9 q = A9 d -A9 u = -1.4Q, 



Qa 



Vu 



0.99746, 



m° t = 150.7 GeV, m° b = 2.5515 GeV, 

1.5142 • 10~ 4 , b u = 0.0395, b' u = 0.0770474, 
q d = 0.0043435, b d = 0.02609, b' d = 0.69138, y d = 0.8100, 
This choice yields at \x = 1 TeV, the following masses and mixings for the quarks: 

1.13 MeV, m c = 0.461 GeV, m t = 150.50 GeV, 
2.53 MeV, m s = 50.99MeV, m b = 2.43 GeV, 



(2.23) 



m d 



\V ( 



CKM\ 



Vw 



0.9745 0.2244 0.0033 \ 
0.2242 0.9737 0.0408 , 
.0.0093 0.0399 0.9991/ 



0.3465, 



(2.24) 



where rjw is the CP violation parameter in the Wolfenstein parametrization. These values, when 
extrapolated to lower energy scales, give extremely good agreement with data [18]. 

We have computed the orthogonal matrices that diagonalize M u and M d . These rotation ma- 
trices will be relevant for our discussion of Higgs-induced flavor violation, as well as FCNC arising 
via SUSY loop diagrams. We find 

/ 0.9645 -0.2640 -0.0001 \ 
-0.1817 0.6642 0.7251 



O u 









0.9840 


-0.1782 


0.0041 ^ 








o d = 




0.1781 


0.9838 


0.0188 










V 


-0.0074 


-0.0178 


0.9998 ) 






I 


0.9988 




-0.0495 


1.17- 10" 








/ 




0.0494 




0.9980 


0.0395 




J 


O u c = 




V 


-0.0020 


-0.0394 


0.9992 


J 






V 



\ 0.1915 -0.6994 0.6886 / 
)988 0.0496 -6.00 ■ 10~ 6 \ 



-0.0494 0.9958 



0.0771 
0.9970 



(2.25) 



/ 



In the case of charged leptons, there is some arbitrariness in the values of (A, B, B', C)e of 
Eq. (11.41) . since we have three observables (charged lepton masses) and four parameters (without 
including the neutrino sector). We shall present a fit with a simplifying assumption B' £ = Be. At 
/i = 1 TeV, a consistent fit for all the lepton masses is found with the following input values: 



A f = 1.67536 GeV, B, = B', = 0.430588 GeV, C e = 0.00742877 GeV 



(2.26) 
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These yield the following eigenvalues at [i — 1 TeV: 

m e = 0.4963 MeV, m M = 104.686 MeV, m T = 1779.5 MeV. (2.27) 

These values correspond to the central values of charged lepton masses when extrapolated down to 
their respective mass scales [18]. The orthogonal matrix that diagonalizes M e is given by 

/ 0.9976 0.0688 9.81 • 10~ 4 \ 
O e = 0.0664 -0.9697 0.2352 , (2.28) 
V -0.0171 0.2346 0.9720 / 

with O e c obtained from the above by flipping the signs in the first row and column. 



3 Symmetry breaking and the Higgs boson spectrum 

We now turn to the discussion of symmetry breaking and the Higgs boson spectrum in the model. 
We shall confine here to the case of having three pairs of Higgs doublets, and no Higgs singlets in 
the low energy theory. It is however, assumed that singlet fields are present in the full theory, so 
that spontaneous Q 6 breaking in the singlet sector appears as soft breaking in the doublet sector. 
As shown in Ref. [8], it is possible to realize such a scenario while preserving the 1 2 interchange 
symmetry for members (1, 2) inside Qq doublets. We seek a consistent picture where CP violating 
phases are generated in the Higgs doublet VEVs. As it turns out, CP also has to be softly broken 
in the bilinear soft SUSY breaking terms, or else there would be no CP phases in the VEVs. 

The superpotential that we consider is the most general one consistent with softly broken Q 6 
symmetry, but preserving the S 2 interchange symmetry: 

W eS = n^Ht + H^H^+n^Hi + ^H^ + H^Hi 

+ pi 31 H»(H« + H d 2 ) + fi 12 {H?H* + H*H%). (3.29) 

As mentioned earlier, we make a 45° rotations in iff, iff and if", ff 2 space, with if" 2 = ^ 2 
and iff 2 = Hl ^ 2 , so that the superpotential becomes 

W eS = (/i! + fi 12 )HiHf + (//1 - ^12)^2 #2 + ^3#3 ^ + V2ii 13 H?H* + V2fi 31 H^ H* . (3.30) 

The redefined fields have {H 2 ) = (H 2 ) = 0. We work in the hatted basis from now on, and drop 
the hat on the new fields. 

The soft SUSY breaking Lagrangian is given, in the rotated basis, as 

U soft = (61 + b 12 )HleHl + (61 - b 12 )H^eH d + b 3 H% eH* 
+V2b 13 H?eH$ + V2b 31 H%eH? + h.c. 

+m 2 dl (\H*\ 2 + \H d 2 \ 2 ) + m%\H d \ 2 + m 2 ul (\H?\ 2 + |if 2 "| 2 ) + m 2 u ,\H%\ 2 , (3.31) 
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where e = ia 2 . 

The full scalar potential including the soft terms, the F terms and the D terms has the form 

V = M 2 dl (\Hf\ 2 + \Hf\ 2 ) + M 2 d3 (\Hf\ 2 + \Hf\ 2 ) 

+ M^(\Hf\ 2 + \Hf\ 2 ) + M 2 3 (|# 3 u0 | 2 + \Hf\ 2 ) 
+ {M 2 3 (Hf*Hf + Hf*Hf) + Ml x (Hf*Hf + Hf* H™ + ) + h.c.} 
+ M 2 d2 (\Hf\ 2 + \Hf\ 2 ) + M 2 u2 (\Hf\ 2 + |tf«+| 2 ) 
+{b[(H? + Hf - HfHf) + b 3 {HfHf - HfHf) 
+ V2b 13 (HfHf - HfHf) + V2b 31 (H% + Hf - HfHf) 
+ b' 2 (H% + Hf - HfHf) + h.c.} 
+ \ g l + gl)(\Hf\ 2 + \Hf\ 2 + \Hf\ 2 + \Hf\ 2 + \Hf\ 2 + \Hf\ 2 
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d -\2 I zj(2°|2^2 



- \ H l \ - \ H 1 I - 1-^3 I ~ I -"3 I ~ \ H 2 I ~ \ H : 

a 2 \m + h 

2 



+-g 2 2 \HfHf + HfHf* + HfHf* + HfHf* 



+H^ + Hf + HfH 2 i *\ 2 . (3.32) 

Here we have redefined new effective parameters for convenience as 

Ml = | /fil+/ x 12 | 2 + 2|/i 3 ir + m 2 1 , M 2 3 = |/i3| 2 + 2|/i 13 | 2 + m 2 3, 

M 2 ul = + /i 12 | 2 + 2|/i 13 | 2 + m 2 ul) M 2 3 = |/i 3 | 2 + 2|/i 3 i| 2 + m 2 3 , 

M d2 = 1^1 - + ^2 = 1^1 - Vl2? + ^ 2 1, 

M 2 3 = y/2(/j,i + Ati 2 )>i 3 + V2/i 3 n* 31 , M 31 = \/2(/xi + ^12)^31 + 

6i = 61 + 612, &' 2 = 61 " 612. (3.33) 

Before analyzing the spectrum, let us note that the potential should be bounded from below 
along all D-flat directions. The following conditions should be satisfied: 

M 2 dl + M 2 ul -2\b' 1 \>0, M 2 1 + M 2 2 >0, M 2 dl + M 2 u3 - 2V2\b 31 \ > 

M 2 d2 + M 2 ul >^ M 2 d2 + M 2 2 -2\b' 2 \>0, M 2 2 + M 2 3 >0, 

M 2 3 + M 2 1 -2v / 2|&i 3 |>0, M 2 3 + M 2 2 >0, M 2 3 + M 2 3 - 2|6 3 | > 0. (3.34) 

In our numerical analysis, we shall verify that these conditions are indeed met. 
We parametrize the VEVs of the four neutral Higgs fields as 

v u i — v sin (3 sin 7„ e tdul , v u3 — v sin (3 cos 7„ e l6uA , 

Vdi = v cos (3 sin 7^ e l6dl , Vds = v cos /3 cos 7^ e t0d3 . (3.35) 



10 



Thus we have |f„i| 2 + | ^3 1 2 + \v d i\ 2 + \vd3\ 2 = v 2 = (174 GeV) 2 . 7^) reflect the orientation of the 
VEVs in the H u r d \\ — H u ^ 3 space, while tan/? is analogous to the up/down VEV ratio of MSSM. 

We can rewrite the potential of the Hi — H 3 sector of the neutral Higgs fields which acquire 
VEVs in a compact form: 



V, 



(1-3) 



N 



HI 



uO* 



(hi 



+ 



+\{gl + gl) 




(10* 



TTdO* 



M 2 dl 

mi: 



M 2 3 
M 2 d3 



dO 



Hf 



+ h.c. 



TTdO* 



do 



Hi 
Hf 



-1 2 



(3.36) 



This suggests a unitary transformation that would diagonalize the first two matrices in Eq. f)3.36p . 
while leaving the D-term unaffected. With such a rotation we have 



V, 



(1-3) 



N 




h% h% 



+ h.c. 



h* h* 
' i 3 'M 



rrin 









(3.37) 



The unitary transformations to go from Eq. (13.361) to Eq. (I3.37P are defined as 





D i{4>u+8 'm-J 



with 9 



A/31 









) 


(3) 




'Hi\ 


)( 





(3.38) 



itan 



-1 



2\Ml\ 



The two phases 0„ and 



M 2 3 - Ml, 
d here are arbitrary. 



- tan 
2 



-1 



2\M 2 3 \ 



M 2 d3 - Ml 



(3.39) 



does not appear in the potential (being 



proportional to U{l)y charges). <p u + <p d can be used to remove one phase of the bilinear terms 
in the potential. Q u ,d are arbitrary diagonal phase matrices. If desired, one can take advantage 
of these phases to remove all but one phase from the parameters of the potential. Since we are 
interested in going back to the original basis from this rotated basis, we find it convenient to set 
Qu,d to be identity. 
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The other parameters of this transformation are 



m 



1,2 



M 2 U , + Mi ± J(M* 3 - Ml,Y + 4IMIJ 



3,1 



M 2 3 + M 2 X ± J(Mj 3 - MIT + A\MU 2 



(3.40) 



and 













= u* v [ 




m 2 24 J 





-b[ -V2b 13 
-V2&31 -63 



(3.41; 



If we choose 



">u + 4>d = 7T + arg 



b'i sino; u sina;rf + V^&si cos u u sin ti^e" 



13 sin uj u cos cj^e 



(3.42) 



m 2 ^ is real and positive (with Q u ^d set to identity). We shall adopt this phase convention in our 
numerical study. However, we shall present analytical results that hold in an arbitrary phase 
convention. 

The task at hand is somewhat simplified, since Eq. (13.371) is relatively simple to analyze. The 



eight real neutral Higgs bosons in can be conveniently parametrized as 

hi = < 

v 3 + 



A5\ 



Vi + -4= (0i + ie(p 5 + ia<fi 7 + i—G) 
V 2 v 



h 3 = e 
^4 



i5s 



+ ibfa + i—G) 
v 



^3 + igfa + icfa -i—G) 
v 



f 4 H — i={4>a + ihcj) 6 + idcj) 7 — i—G). 
V2 w 



(3.43) 



Here Vi (i — 1, 2, 3, 4) are the magnitudes of the VEVs of the redefined fields /ij, and Si are their 
phases. Without loss of generality we have taken t> 4 to be real. G in Eq. (13.431) is the Goldstone 
field eaten up by the Z gauge boson. We shall work in the unitary gauge and set G = 0. We have 
checked explicitly that the G field does not mix with other scalar fields, and that its mass is exactly 
zero. The coefficients of various fields in Eq. (I3.43[) are functions of the Vi's: 



V3 



2 

4 b 



f 



V2y/vi+i 
v^v\ + v. 

V V 2 + V l ' 



^3\Ai + v. 



d _ VW V 1+ V 3 



h 



V2 



(3.44) 
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We shall allow for the soft SUSY breaking parameters (pi) in the Higgs potential to be complex. 
Phase rotations cannot remove all phases from the potential, one phase is unremovable. Without 
this phase, the model cannot induce complex VEVs to the doublets, as shown in Ref. [19] by 
a geometric argument. For the case when all parameters in the Higgs potential are real, we have 
numerically verified that the CP violating extremum would generate two massless modes, signalling 
inconsistency with symmetry breaking [19]. 

We take the soft bilinear terms mf 3 , m^ 4 , m 23l m 2A of Eq. (I3.37P to be complex, and denote 
the phase of mf„- as The minimization conditions then read as 

m\vi + \m 2 13 \v 3 cos{6i 3 + 6 1 + S 3 ) + \m 2 u \v 4 cos(8 u + St) + -^{gl + sfMC^i + v l ~ v l ~ v l) = °> 

m\v 2 + \ml 3 \v 3 cos(8 2 3 + S 2 + S 3 ) + \ml^\v 4 cos(<5 2 + 9 24 ) + ~^g\ + g 2 ) v 2{v\ + v\-vl- v\) = 0, 

m 2 3 v 3 + \ml 3 \vi cos(6»i 3 + 6 1 + 8 3 ) + \m 2 23 \v 2 cos{6 23 + 5 2 + 8 3 ) - i(</J + g 2 )v 3 {vf + v 2 - v 2 3 - v\) = 0, 

mlv4 + \ml 4 \vi cos(6 l i4 + St) + \ml^\v 2 cos(5 2 + #24) - ^(#1 + #2)^4(^1 + v\ - v\ - vl) = 0, 
\m\ 3 \{v\ + vl) sin(6 ) i 3 + Si + S 3 ) + \m 23 \viv 2 sin(6 , 2 3 + 5 2 + S 3 ) + \m 2 4 \v 3 V4sm(9 14 , + Si) = 0, 
\ m 2A\( v 2 + v l) sin(6> 24 + S 2 ) + \mu\viv 2 sm(9 u + St) + \ml 3 \v 3 v 4 sm(6 23 + S 2 + S 3 ) = 0, 
|mi 4 |fit;4sin(6 l i4 + St) - \m\ 3 \v 2 v 3 sin(9 23 + S 2 + S 3 ) = 0. (3.45) 

Denoting the squared matrix for (pi, i = 1,2,.. .7 from the Hi — H 3 sector as 

^0,(1-3) — J^lj 

we obtain 



Kd-3) =Ml, (3.46) 



M 2 n = Xv l + K^[cot(^ 4 + St) - COt (5i3 + Si + S 3 )], 



V 2 



M 2 22 = \V 2 + K [cot (^23 + S 2 + S 3 ) - COt(# 24 + S 2 )), 
V 2 



M 2 33 = Xvl + «^[cot(0 23 + S 2 + S 3 ) - COt(^i 3 + St + S 3 )), 

V 3 

M\ A = Xvj + K— [cot (014 + St) - COt(#24 + S 2 )), 
l>4 

M ^5 = K ^TT^ft COt (^14 + Sl) + V ~\ COt (023 + S 2 + 5 3 ) - ^ t 2^ COt (^13 + <*1 + $0], 
Vf + W3 W{ fg V{V 3 

M W = K -^l V l COt (^14 + Si) + V\ COt(0 23 + S 2 + S 3 ) - {v\ + Vl) 2 COt(0 2 4 + S 2 )}, 

V 2 V4 : {V 2 + V 4 ) 



= K V2V f^+il\ ; j + ^ 2 ) [cot(g + 5i) + cot( ^ + 52 + 5a)] 
K +^3)^2 + ^4) " 



A<i2 = Af 1^2, A4? 3 = -Awiw 3 + K,——— cot(6'i 3 + <5i + 5 3 ), 



13 



M 2 



M 



14 



2 

17 



M 2 



24 



A* 



2 

27 



M 2 



36 



A* 



2 
46 



M 2 



Ml r = 



M 2 



67 



V 2 

-\V1V4 — K — COt (^14 + S\), 
Vl 



M 2 15 = 



-K- 



V2V4 



V 3 



Mj 6 = K 



v 2 



K- 



V 2 V A 



vS 2 + u 



Vl y/vf + vly/vfT 



M 



2 

23 



—\v 2 V 3 - K COt (0 2 3 + S 2 + S3), 
^3 



-\v 2 V A + KCOt(0 2 4 + S 2 ), 
V 4 



M 2 



25 



K- 



V\V A 



V 3 



12 j 



M 



2 

26 



— K 



V 2 



T2' 



— K- 



«2 



«2 



W 4 



Al 2 



37 



— K 



Af 34 = A-y 3 t> 4 , Al 
^2^4 y/vl + vl + vl + v. 



2 

35 



— K- 



V 2 V 4 



— K 



1 



M 



2 

47 



^3 y/vj + V%y/v% + v\ 

v 2 \/v\ + 



Al 



T2' 



^3 



vi \Jv{ + u| 



*2 



= K 



K 



= K 



V1V3 \/v\ + vl yjv\ + V\ 
V 2 V 4 \/v\ + v\ + V 2 + V, 
VlV 3 {vl + vl)y/vl + V\ 

^v\ + vl + vl + v. 



V v f + v iV v 2+vl 

v\v\ COt (014 + 81)+ v\v\ COt (023 + S 2 + S3)], 



■[vl COt (014 + Si) - v\ COt (023 + S 2 + <J 3 )], 



2 

4 L.2 



v{ cot (0i4 + <5i) - U4 cot(0 23 + 5 2 + 5 3 )]. 



V^F + ^I(^2 2 + vf) 

Here we have defined A = (gj + <?|)/2 = M 2 z jv 2 and k = m| 4 sin(0 2 4 + #2)- 
The potential of the H 2 — H 2 fields which do not acquire VEVs is 



(3.47) 



V, 



(2) 



N 



M 2 2 \H%\ 2 + M 2 2 \H d 2 \ 2 - {b 2 H-H d 2 + h.c.} 



+ 



gl + gl 



l\ ttu 1 2 I Tjd\2 , I 12 1 I 12 I 12 
(|i7 2 | -|i? 2 | +Kl| + |U«3| - |Vdl| -\Vd3 



2\2 



(3.48) 



The corresponding squared mass matrix for the scalars in the basis (ReH 2 , lmH 2 , ReH 2 , lmH 2 ) 



is 



/ M 2 2 -^cos2/3 



M 2 0{2) 



\ 



2 


Reb' 2 
—hnb'n 



M 2 2 ~^cos2(3 
—lmb' 2 
—Reb'r, 



Reb 2 
— Imb'o 



M j 2 + ^ CO s2(3 



— Im& 2 
— Reb' 2 


M 2 2 + ^cos2/3 J 



. (3.49) 



This matrix has two pairs of degenerate eigenstates, owing to an unbroken U(l) symmetry. 

The H^ ,d —H^ ,d sector charged Higgs boson mass matrix is, in the basis {H™ + , H% + , Hf , H$ }, 



with 



M 2 U = M 2 



ul 



^1(1-3) 



-m| cos 2(3 + hgl \v dx \ 2 , 



(M 2 )ij, 



M 2 22 = M, 



u3 



-m 2 z cos 2(3 + -gj\v d3 \ 2 , 
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M 2 33 = M 2 dl + im| cos 2/3 + ^g 2 2 \v ul \ 2 , M 2 U = M 2 3 + cos 2/3 + ^g 2 2 \v u3 \ 2 , 
M\ 2 = M\* = M 2 31 + M\ 3 = M\* = b[ + ^g 2 2 v* ul v* dl , 

M 2 U = M 2 41 * = v^&is + 2& v >di> M lz = ^'2* = V2&31 + ^Xi^ 

M 2 24 = Ml* = b 3 + ~s 2 2 < 3 <4, M 2 M = Ml,* = M 2 3 + \g\v ul vl 3 . (3.50) 

Finally, the — H 2 sec t° r charged Higgs mass matrix is, in the basis {-ff 2 + , }, 



U 2 I M 2 2 - ±m|cos2/3 b 2 , 

M ^-\ b ' 2 * M 2 2 + Im|cos2/3 1 



Now we present two sets of numerical fits (cases (1) and (2)) which show the consistency of 
symmetry breaking. We are interested in choosing the SUSY breaking parameters (including the 
H terms) around the TeV scale, guided by arguments of naturalness. At the same time we wish 
the spectrum to be consistent with FCNC constraints arising from meson-antimeson mixings. We 
have explored parameter space of the Higgs potential where both these constraints are met. For 
the FCNC constraint, we allow the new Higgs exchange contribution to AM be not more than the 
experimentally measured values. 
Case (1) 

The parameters in the original Higgs potential of Eq. (I3.36P are taken to have the following 
values. 

M dl = 3.754 TeV, M d3 = 3.586 TeV, M ul = 4.782 TeV, M u3 = 2.152 TeV, 

M 31 = 2.336 e i0 - 792 TeV, M 13 = 1.346 e - iL205 TeV, b[ = 3.144 e i2 - 963 TeV 2 , 

b 3 = 3.196 e i2 - 064 TeV 2 , b 31 = 4.052 e i2A86 TeV 2 , b l3 = 3.438 e i3 - 109 TeV 2 , 

M u2 = 4.550 TeV, M d2 = 4.850 TeV b' 2 = 0.000 TeV 2 , (3.52) 

In the representation of Eq. (13.37|) this choice corresponds to 

mi = 4.937 TeV, m 2 = 1.767 TeV, m 3 = 3.923 TeV, m 4 = 3.401 TeV, 

m 13 = 1.851e- iL437 TeV, m u = 2.736 e - ia732 TeV, m 23 = 2.442 e iL347 TeV, 

m 24 = 2.104 TeV . (3.53) 

For completeness we also give values of other parameters, u u = 0.70, u d = 0.622, (p u + (p d = 1.005. 
We obtain numerically the VEV parameters to be 

tan/3 = 2.00, A8 d = -0.03, A6 U = 1.37, tan 7d = 2.50, tan 7u = 0.33. (3.54) 
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The mass eigenvalues of the Higgs bosons in the Hi — H 3 sector are found to be 

M h0 = (99.4, 115.1) GeV, M x = 3.299 TeV, M 2 -M x = 0.226 GeV, 
M 3 = 4.161 TeV, M 4 — M 3 = 0.411 GeV, M 5 = 5.124 TeV, M 6 - M 5 = 0.040 GeV. (3.55) 

Note the appearance of nearly degenerate states (Mi, M 2 ) etc, with their mass splitting being 
proportional to m|/4. The Higgs bosons from the H 2 sector have degenerate masses given by 

M 7 = M 8 = 4.850 TeV M 9 = M 10 = 4.550 TeV. (3.56) 

The charged Higgs bosons are nearly degenerate with its neutral partner, so we list the mass 
splittings: 

Af ± i - Mx = -0.532 GeV, M± 2 — M 3 = —0.156 GeV, M± 3 — M 5 = 0.032 GeV . (3.57) 

In the [H 2 — H 2 ) sector, the two charged Higgs bosons are degenerate with the neutral ones 
given in Eq. (13. 56ft . 

The mass eigenstates Hi are mixtures of hi, i — 1, 2, ... 7 states in the (1-3) sector. The 
orthogonal transformation that diagonalizes the mass matrix of Eq. (13.461) is 



/ 


0.0662 


0.8919 


0.2708 


0.3562 


8.60 ■ 10~ 7 


1.23 ■ 10~ 6 


2.15 ■ 10" 6 \ 




f 


h 


\ 




0.0314 


-0.0023 


0.0427 


-0.0324 


-0.4002 


0.8800 


0.2482 






ho 






0.3322 


-0.2620 


-0.3269 


0.8428 


0.0204 


0.0293 


0.0509 






h 






-0.0357 


0.0026 


-0.0484 


0.0368 


0.1514 


0.3354 


-0.9272 






hi 






-0.0644 


0.3645 


-0.9010 


-0.2159 


0.0231 


0.0332 


0.0578 






h 5 






0.0430 


-0.0032 


0.0584 


-0.0444 


0.9029 


0.3311 


0.2607 






h 6 




\ 


-0.9365 


-0.0553 


-0.0289 


0.3345 


0.0279 


0.0401 


0.0697 ) 






h- 


J 



(3.58) 

with k = 0, • • ■ 6. Since b' 2 = in this case, the H 2 mass matrix is diagonal, and thus the mass 
eigenstates are the original state. 
Case (2) 

Here we take the input parameters corresponding to Eq. (13.36!) to be 

M dl = 3.980 TeV, M d3 = 5.412 TeV, M ul = 2.765 TeV, M u3 = 3.692 TeV, 

M 31 = 2.825 e ia781 TeV, M 13 = 1.693 e " ia949 TeV, b[ = 3.698 e ilA95 TeV 2 , 

b 3 = 3.097 e iL522 TeV 2 , 6 3 i = 7.420 e i2A28 TeV 2 , 6i 3 = 1.840 e " i2 772 TeV 2 , 

M u2 = 3.550 TeV, M d2 = 5.850 TeV, b' 2 = 1.234e iL56 TeV 2 . (3.59) 
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This choice corresponds to parameters in Eq. (13.371) to be 



mi = 4.377 TeV, m 2 = 1.154 TeV, m 3 = 5.466 TeV, m 4 = 3.906 TeV, 
m 13 = 3.281e iL2n TeV, m u = 1.702e ia974 TeV, m 23 = 3.190 e ~ ia501 TeV, 
m 24 = 2.326 TeV, 



(3.60) 



with u u = —0.501, uJd 



-0.606, (p u + (pd = 4.786. 



The Higgs VEV parameters are found for this input to be 

tan/3 = 2.40, A6 d = -0.06, A9 U = 1.34, tan j d = 1.80, tan 7u = 1.00. 
The mass spectrum of Higgs boson in the Hi — H 3 sector is 



(3.61) 



M h0 = (104.1, 119.2) GeV, M l = 2.869 TeV, M 2 - M l = 0.325 GeV 
M 3 = 5.114 TeV, M 4 - M 3 = 0.132 GeV, M 5 = 5.658 TeV, M 6 - M 5 = 0.087 GeV, (3.62) 
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3.541 TeV. 



while the mass eigenvalues of Eq. (13.491) are 

M 7 = M 8 = 5.856 TeV M 9 = .U 
The charged Higgs boson masses are given by 

M±i — Mi = 0.225 GeV, M± 2 — M 3 = 0.182 GeV, M ±3 - M 5 



-0.064 GeV . 



(3.63) 



(3.64) 



with the remaining two charged Higgs bosons being degenerate with the neutral ones given in Eq. 

dSSSD- 

The orthogonal matrix that diagonalizes Eq. ( 13.461) is 



H 



/ 


0.3919 


0.8356 


0.2620 


0.2819 


1.02 • 10~ 4 


-5.99- 10~ 5 


7.35 • 10~ 5 


\ 


/ 


hi 


\ 




0.4761 


-0.2234 


-0.1898 


0.1764 


0.6693 


0.4070 


0.2065 






h 2 






-0.5233 


-0.0380 


0.4553 


0.4166 


0.4223 


-0.2386 


0.3295 






h 






0.5065 


-0.2376 


0.1014 


-0.0942 


0.0903 


-0.8105 


0.0535 






hi 






-0.0893 


0.2828 


0.0490 


-0.7599 


0.1969 


0.0134 


0.5415 






h 5 






-0.0625 


0.0293 


0.3803 


-0.3534 


0.4512 


0.0437 


-0.7212 






h e 




\ 


0.2783 


-0.3363 


0.7285 


-0.0671 


-0.3510 


0.3422 


0.1805 


/ 


\ 


h- 


J 



(3.65) 



Eq. (|3.49| 


is 




-0.0568 


0.000 


0.9984 


0.0000 - 


-0.0568 


0.0000 


0.0000 


0.9984 


0.0000 


0.9984 


0.0000 


0.0568 



V Im(^ 2 d ) / 



(3.66) 
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with k = 7, • • • 10. 

In these fits, Mh is the light standard model-like Higgs boson mass, for which radiative correc- 
tions are significant. In our computation we have included known two loop corrections. The two 
values listed for Mh correspond to zero and maximal left-right stop mixing (X t = or 6). We 
have taken m t = 174 GeV, Msusy = 1.5 TeV and a s (m t ) = 0.108 for these evaluations and used 
the analytic approximation given in Ref. [20]. 

An interesting feature of these two fits is that the diagonal entries of the quadratic mass matrix 
of the potential of Eq. (I3.36f) are all positive. This of course does not preclude some soft squared 
masses turning negative as in the MSSM via large top quark Yukawa coupling (since the diagonal 
entries also receive /i term contributions), however, this is not necessary for symmetry breaking 
to be triggered. Yet, one of the eigenvalues of this matrix is negative, which facilitates symmetry 
breaking. For the two cases we find these eigenvalues to be 

Case (1) : {(5.123 TeV) 2 , (4.161 TeV) 2 , (3.300 TeV) 2 , - (38.682 GeV) 2 } , 

Case (2) : {(5.658 TeV) 2 , (5.115 TeV) 2 , (2.869 TeV) 2 , - (45.40 GeV) 2 } . (3.67) 

The conditions for boundedness of the potential listed in Eq. ( I3.34p are found to be satisfied for 
both cases. 



3.1 Neutralino and Chargino masses 



The symmetry breaking parameters do not fully determine the masses of the neutralinos and the 
charginos. Here we present analytical results for their mass matrices 
The mass matrix of Hi 



~ ~ n ~o ~o ~o ~o 
H 3 sector neutralino in the basis of {B, W v , H% , H% , Hf , H$ } 



is 



M v o 



X°(13) 



The mass matrix of the H 2 sector in the basis { H% , H$ } is 







X°(2) 




/ 


M 6 





92 f ul 


92fu3 


92V dl 


92V ds 


\ 


V2 


V2 


V2 


V2 







M W 


9i fui 


9lf«3 


gifdi 


9lf<i3 






V2 


V2 




V2 






92 Vul 
V2 


91 fui 








-(fix + H12) 








92V 
V2 


giv U 3 

V2 








-V2fi 3 x 


-Ms 






92V dl 
V2 


givdi 

V2 


— (fll + lilt) 












V 


92V dZ 
V2 


9lVd3 
V2 




-/is 








/ 



(3.68) 



(3.69) 



The mass matrix of charginos of the Hx—H 3 sector in the basis {IV + , H\ , H$ , W~ , Hf , H% } 
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has a block-diagonal form: 



M 



X±(13) 



X T 
X 



with 



X 



/ M w g x v ul g x v u2> ^ 
\ giv d3 V2fi 13 fi 3 J 



The chargino mass matrix in the H% — sector in the basis of , } is 

M x ± {2 ) = 



/Ii - /ii2 

Mi - Hl2 



(3.70) 



(3.71) 



(3.72) 



4 Tree level Higgs induced FCNC processes 

In this section we discuss various FCNC processes mediated by tree-level neutral Higgs boson 
exchange. 



4.1 Neutral meson mixing via Higgs exchange 

Accurate measurements exist [21] for neutral meson-antimeson mixings in the K° — K°, B® — B®, 
B® — B® and in D° — D° sectors. In the Qq model there are new contributions to these mixings 
arising through tree-level Higgs exchange. These new contributions will modify the SM predictions, 
which are all in good agreement with data. Here we compute these new contributions, following 
the analysis of Ref. [22], with updated QCD corrections and hadronic matrix elements. 

The Yukawa coupling a Ui f3 u , a, U & *V d °f Eq. d2HlJ) can be determined from the mass 
matrix Eq. f[2~lH : 



\Vu, ctt\ 



m 



t, b u, d 



\V U , dl\ 



y, d 



_ ml b b U) d 

d — — j— j— 

l^u, rfl | 

m l bVu, d/y u , d 

\ v u, ds\ 



(4.73) 



Using the input values given in Eq. (12.231) we get for the two cases 
Case (1) 



Old 



0.0409, 
= 0.7195, 



Pd~- 

Pu 



6.51 • 10" 4 , p' d = 0.0173, 5 d = 3.35 • 10~ 4 , 



0.0858, 



X = 0.1672, 5 U = 1.10- 10~ 4 . 



Case (2) 



19 



a d = 0.0526, (3 d = 7.46 ■ 10" 4 , d = 0.0198, 5 d = 4.30 • 10~ 4 , 
a u = 0.9354, (3 U = 0.0372, U = 0.0724, 5 U = 1.43 ■ 10~ 4 . 
After 45° rotation in the Q§ doublet space, the Yukawa coupling matrices in the down sector 



arc 



Y, 



di 



T d P d 



( \ 

(3 d 
\ d J 



Y, 



,12 



T d P d 



I o o p d \ 


\0 d J 



Pd,cO d c, 



Y 



d3 



T d P d 



Pd c O d c, 



(4.74) 



/ 5 d \ 

-s d 

\ a d / 

where P d , P d c are defined in Eq. (12.181) . The Yukawa couplings in the up-quark sector and the 
charged lepton sector are similar. 

The new Higgs-mediated contributions to AF = 2 Hamiltonian, responsible for the neutral 
meson-antimeson mixings has the form [22] 

2 



H. 



1 



cff 



2M h 2 



yk 1 + 75 y k * 1 - 75 
ij O J* 9 



9j 



(4.75) 



Here qij are the relevant quark fields contained in the meson. are the Yukawa couplings of qi, 
qj with Higgs mass eigenstate H k mediating FCNC interactions, k — 1, 2, ... 10 in our model, 6 
from the (Hi — H%) sector and 4 from the H 2 sector. (The light standard model-like Higgs boson 
has practically no FCNC couplings.) Y^ can be obtained via inverse transformations, Eq. ( 13. 38ft . 
(EH and (GLSED or (l3~65l) . 
We obtain 

,2 



Mf 2 = (<P\H eS \ 



2M, 2 



24 (m qi + m q .) 2 



4-Y k Y k * 

' ij ji 



y fc2 _i_ y k * 2 

ij 



12 2(m qi + m qj ] 



B 2 -7)2(0) 
■ B 4 ■ T] 4 (fJ,) 



(4.76) 



Here <fi is the neutral meson (K°, B^, B° s , D°). For our numerical study we use the modified 
vacuum saturation and factorization approximation results for the matrix elements [2, 3] 

(0|/,(1± ls)fji(l Tls)fj\ 



V6 (m qi + m qj y 



■ B 



(Pf,(\ : '.-,)/,./;( 1 : '.-J/; 



5,2 ^ 

6 (m qi + m qj y 



B, 



(4.77) 



i?2 and .84 are equal to one in the vacuum saturation approximation, but are found to be slightly 
different from one in lattice simulations. We use (B2, B4) = (0.66, 1.03) for the K° system, 
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(0.82, 1.16) for the Bj and B° s systems, and (0.82, 1.08) for the D° system [2]. In Eq. ffl~T6l) rjM, 
rj^(fi) are QCD correction factors of the Wilson coefficients C 2 and C4 of the effective AF = 2 
Hamiltonian in going from the SUSY scale M s to the hadronic scale \x. These factors are computed 
as follows. The AF = 2 effective Hamiltonian has the general form 

5 3 

ntl =2 = ^Ci Qi + ^Ci Q, (4.78) 

i=l i=l 

where 

Qi = qili^j a L qih U( ijL^ Q2 = ftS^**^ Qs = aS?*??**^ 

Q4 = qi R qj a L qdq jR , Q?> = ftS<&f«J<&£, (4.79) 

with Qi,2,3 obtained from Qi.2,3 by the interchange L <-> R. 

For computing 772,4 we take the SUSY scale M s to be 1 TeV. All the supersymmetric particles 
and heavy Higgs bosons are integrated out at 1 TeV. The Wilson coefficients evolve from M s down 
to the hadron scale \i according to the equations 

CM = E J2^' S) + V^ s) )v a *C s (M s ), (4.80) 

i s 

Here r] is defined as rj = a s (M s )/a s (m t ). The magic numbers Oj, bf"'^ and c[ ,s ^ can be found in 
Ref. [2] for the K system, in Ref. [3] for the Bd >s system and in Ref. [4] for the D system. With 
M s = 1 TeV and a s {m z ) = 0.118, and m t {m t ) = 163.6 GeV we find r\ = a s (l TeV) / a s (m t ) = 
0.0882/0.108 = 0.8167. 

At the SUSY scale, the neutral Higgs bosons in our model generate only operators Q2 and Q4. 
Consequently, at the hadron scale, for the K° system, we find 

C 2 (fi) = C 2 (M S ) ■ (2.54), CM = C 4 (M S ) ■ (4.81), 
CM = C 2 (M S ) ■ (-1.8 x 10- 3 ), CM = C 4 (M S ) • (0.186), (4.81) 

leading to r] 2 {^) = 2.54, ^(/x) = 4.81. Although operator mixings induce non-zero C3 and C5 at 
the hadronic scale, their coefficients are found to be rather small. 
For the B® s system, following the same procedure, we find 

CM = C 2 (M S ) ■ (2.00), CM = C,{M S ) ■ (3.12), 
CM = C 2(M S ) ■ (-2.44 x 10" 2 ), CM = C4(M S ) • (0.0874). (4.82) 

And for the D° system we have 

CM = C 2 (M S ) ■ (2.31), CM) = C A {M S ) ■ (3.99), 
CM = C 2 (M S ) ■ (-1.30 x 10- 2 ), CM = C4(M S ) • (0.144). (4.83) 
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In all cases we see that the induced operators C 3 and C 5 are negligible. 

K° — K° mixing constraint: 

In the K° system, tree-level neutral Higgs boson exchange contributes to Kl — K$ mass 
difference, as well as to the indirect CP violation parameter, modifying the successful SM pre- 
dictions. The mass difference is computed from Aitlk = 2ReM^, while the CP violation pa- 
rameter is | | — ^=^-. We seek consistency with the precisely measured experimental val- 
ues Amx/iTLx — (7.1 ± 0.014) x 10 -15 and |e^| ~ 2.3 x 10~ 3 . In our calculation, we choose 
rriK = 498 MeV and fx = 160 MeV. For the two numerical fits we find the new contributions to 
be 

Case (1) : (Am K /m K ) ncw = 7.361 x 10~ 15 , e^ cw = 2.00 x 10~ 4 , 

Case (2) : (Am K /m K ) ncw = 5.721 x 10~ 15 , e^ cw = 2.28 x 10~ 5 . (4.84) 

The contributions from H® — sector and H® sector to Re(M^) are respectively (3.033 x 10~ 15 , 
-1.200 x 10~ 15 ) GeV for case (1) and (2.512 x 10" 15 , -1.088 x 10~ 15 ) GeV for case (2). We see 
that the new contributions to the mass difference is significant, but consistent with data. New 
contributions to CP violation is suppressed, which is a generic feature of Higgs exchange in this 
class of models. We elaborate on this issue later in this section. 

B® — mixing constraint: 

For the B®—B° d system We use as input m Bd = 5.281 GeV, fs d = 240 MeV and seek consistency 
with the experimental value Am Bi = 3.12 x 10~ 13 GeV. We find for the Higgs induced contribution 

Case (1) : (Am B J ncw = 2.997 x 10~ 13 GeV, 

Case (2) : (Am Bd ) ncw = 2.728 x 10~ 13 GeV. (4.85) 

The contributions from H%-H% sector and H% sector to M h £ are (2.298 x 10~ 14 , 1.269 x 10~ 13 ) GeV 
for case (1) and (2.137 x 10~ 14 , 1.150 x 10~ 13 ) GeV. Again, we see consistency with experimental 
values. CP violation parameter is found to be extremely tiny, ~ 10~ 5 , from Higgs boson exchange. 

B° s — B® mixing constraint: 

For the B° s — B® system, we use rag, = 5.37 GeV, f Bs = 295 MeV and compare the new 
contributions with Am Bs = 1.067 x 10~ n GeV. 

Case (1) : (Am B J new = 1.688 x 10" 12 GeV, 

Case (2) : (Am B J ncw = 1.396 x 10~ 12 GeV, (4.86) 

The H^-H$ sector and the H$ sector contribute to M§ given by (8.532 x 10~ 13 , -9.460 x 10~ 15 ) 
GeV for case (1) and (7.067 x 10~ 13 , -3.835 x 10" 15 ) GeV for case (2). These new contributions 
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are within experimentally allowed range. Higgs mediated CP violation is again found to be highly 
suppressed. 

D° — D° mixing constraint: 

For the D° — D° mixing we use mo = 1.864 GeV, fo = 200 MeV and compare the new 
contribution with Am D = 1.27 x 10 -12 GeV. 

Case (1) : (Am D ) new = 8.620 x 10" 13 GeV, 

Case (2) : (Am D ) ncw = 2.645 x 10" 13 GeV, (4.87) 

The H® — H® sector contribution has different sign from that of the H® sector. We find for 
these contributions to be (4.402 x 10~ 15 , -4.354 x 10~ 13 ) GeV for case (1) and (2.568 x 10~ 15 , 
— 1.348 x 10~ 13 ) GeV for case (2). Again these limits are within experimental range. 

We have found that new sources of CP violation through tree-level Higgs is very small in 
meson-antimeson mixings with typical values Im(Mi2) ~ 10~ 4 Re(Mi2). This can be understood 
heuristically as follows. There are two types of contributions to the meson mixing as given in Eq. 
(14.761) . The first term, proportional to B 2 respects a global U(l) symmetry (strangeness in the K° 
system) , which is only broken by the mass-splittings in the neutral Higgs boson spectrum between 
a pair of particles. However, this splitting is very small, of order m\ in the squared mass, see Eqs. 
(13.551) . The couplings of the nearly degenerate Higgs in each pair differ by a factor i, owing to 
the U(l) symmetry, and the two contributions cancel, in the limit of exact degeneracy. For both 
the real and imaginary parts of M 12 the contributions from the first term is suppressed by a factor 
m|/(4M|). Such a suppression is absent in the second term of Eq. (14.761) . since the operator Q4 
explicitly breaks the U(l) symmetry. Thus, although the first term has CP violation, in relation to 
the CP-conserving second term, it is suppressed by a factor m 2 z / (AMf:) ~ 10~ 4 . Now, the second 
term, while it has no suppression factor, it is purely real. This can be seen from the following 
observation. In the mass basis of fermions in the original basis we have the relation (owing to the 
vanishing of off-diagonal mass terms in the mass eigenbasis) 

(Y^ (H$) = (H*) (4-88) 

for i 7^ j. The couplings of mass eigenstates of Higgs boosn to down-type quarks are simply linear 
combinations of Hf and . Since we assume CP to be spontaneously broken, all components of 
(Y k )ij with i 7^ j have the same phase. As a result the second term of Eq. 14.761 becomes real. The 
constraint imposed by SUSY, that H* fields do not couple to down-type quarks, and the fact that 
only two of the down-type Higgs bosons acquire VEVs is very crucial for this result. 
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4.2 Neutron electric diploe moment from Higgs exchange 

Higgs boson exchange can generate non-zero electric dipole moments for the fermions. These 
diagrams are however suppressed by the light fermion Yukawa couplings. For the d quark EDM 
arising from neutral Higgs boson exchange at the one-loop level we find [22] 



-In 



Ml 
m 2 , 



id. 



(4.89) 



where ^ = (a s (M k ) /a s (fx)) 16 ^ 23 0.12, and q is summed over d, s and b. The neutron EDM is 
determined using the quark model via 



D n = Ad d /3 - d u /3. 



(4.90) 



We find 



Case (1) : D n = 1.809 x 10 



-31 



cm, 



Case (2) : D n = 6.091 x 10~ 31 e - cm, 



(4.91) 



which are well within experimental limits. The EDM of the electron is similarly found to be 
extremely small from the Higgs boson exchange diagrams. 



4.3 fi — > 3e and r — ► 3/i decays 

Tree-level Higgs boson exchange can lead to flavor violating leptonic decays such as r 
fx — > 3e. The effective weak interaction mediating such decays can be parametrized as 



G, 



eflf 



E(^)n(n) 



Mi 

3e for the two cases to be 



The effective couplings are found for fx - 

Case (1) : G cS = 4.432 x 10~ 13 G F , 
Case (2) : G cS = 4.191 x 10" 13 G F . 

3 fx decay are 



And the couplings for r 



Case (1) : G cS = 45.721 • 10~ 8 G F , 
Case (2) : G cS = 6.977 • 10" 8 G F . 



3 fx and 



(4.92) 



(4.93) 



(4.94) 



Such small effective couplings will lead to negligible contributions to the decay branching ratios. For 
example, the branching ratio for r — > 3fx is of order 10~ 15 , well below the experimental sensitivity. 
We conclude that Higgs mediated FCNC in the lepton sector are all safe. 
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5 FCNC mediated by SUSY particles 



In this section we turn attention to the flavor changing processes mediated by the supersymmetric 
particles. The main motivation for the non-Abelian Q$ model was to bring such processes under 
control by a symmetry reason. Here we analyze meson-antimeson mixings, flavor violating leptonic 
decays, and the EDM of the neutron and the electron. We present our proposal to suppress SUSY 
contributions to the EDM by making the Higgsinos of the model light, with masses of order 100 
GeV. 

Owing to the Qq symmetry, the first two family squarks (and similarly sleptons) are degenerate 
in mass, while the third family, which is a Q e singlet has a different mass. In the fermion sector 
Q 6 symmetry is broken, which means that there will be SUSY loop induced flavor violation in the 
model. Constraints on such flavor violation has been listed in Ref. [2-4] assuming all three families 
of squarks are degenerate. While these results are applicable for the K° and D° system in our 
model, they do not work well for the B® s system. This is because the masses of the b and d, s 
masses are not the same. 



5.1 Generalized constraints for Bd system 

We have generalized the results of Ref. [3] by allowing for b mass to be different from the masses 
of d, s. We define new parameters 

y% B = — 2 — (5-95) 

m dA,B 

for A, B = L, R. We expect these y parameters to be of order one, but not very close to one. 
Taking account of y ^ 1 we have generalized the constraints on the squark mixing parameters from 
B° d system as follows. 

The effective AF = 2 Hamiltonian for B d:S system can be written as 



'Hcs = ^ Ci Qi + ^ C-i Q 



i=i i=i 

!K {(5 d 13 )l L (2AQ 1 xf e (x,y) + 66Qj e (x,y)) + (5 d 13 ) 2 RR (2AQ 1 xf 6 (x,y) + 66Q 1 Ux,y)) 



216m 2 - 

a 



+(S d 13 ) LL (5 d 13 ) RR (50AQ 4 xf 6 (x,y) - 72Q 4 f e (x,y) + 24Q 5 xf 6 (x,y) + 120Q 5 f 6 (x,y)) 
+(8 d 3 ) 2 RL (20AQ 2 xf 6 (x,y) - 36Q 3 xf 6 (x,y)) + (8 d 3 ) 2 LR (20AQ 2 xf 6 (x,y) - 36Q 3 xf 6 (x,y)) 
HS d 13 )L R (5 d 13 ) 2 RL (-132Qj(x,y) - 180Q 5 f 6 (x,y))}, (5.96) 

The functions f G (x,y) and f&(x,y) are 

f 6 (x, y) = (a ,_ 1)3 ( y ^ 1 )3( ;z _ 1 )3 [ - ln x ( x + y + % y - 3x ^(y - x ) 3 
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+ \ny(x + y + xy - 3y 2 )(x - l) 3 

+2(x -l)(y- l)(-x + y + x 2 - y 2 - x 3 + y 3 + 2x 2 y - 2xy 2 )^ 

/«(*» V) = {x - my - mz -i r In x • x(a; 2 — — l) 3 
+2\nyy(x-y 2 )(x-l) 3 

+(x -l)(y- l)(x 2 - y 2 + x 3 - y 3 - 7x 2 y + 7xy 2 + x 3 y - xy 3 ) . (5.97) 

Generalizing the results of Ref. [3] we obtain the squark mixing coefficients (<^ 3 )ab with A, B = 
(L, R) as shown in Table 1. Here we have used the same input as in Ref. [3], so that for y = 1 our 
results coincide. We have used the next-to-leading order lattice calculation results for the matrix 
elements. For some of the mixing parameters we made a simplifying assumption that yf and y^ 
are equal. 



y\x 


0.25 


1.0 


4.0 


0.25 


1.0 


4.0 




|Re(^ 3 ) LL 




Im(<J? 3 )LL 




0.25 


3.4 x 10" 2 


1.6 x 10" 1 


2.5 x 10" 1 


7.2 x 10" 2 


3.4 x 10" 1 


1.2 x 10" 1 


1.0 


6.2 x 10~ 2 


1.4 x 10" 1 


7.0 x 10- 1 


1.3 x 10" 1 


3.0 x 10" 1 


3.4 x 10" 1 


4.0 


1.6 x 10- 1 


2.7 x 10- 1 




3.3 x 10- 1 


5.8 x 10- 1 






|Re(^ 3 ) RR 


= |Re(5f 3 ) LL 




|Im(^ 3 ) RR 


= |Im(^ 3 )LL 




0.25 


1.4 x 10~ 3 


2.4 x 10" 2 


1.0 x 10" 2 


4.4 x 10" 3 


1.0 x 10~ 2 


4.3 x 10~ 3 


1.0 


1.9 x 10~ 2 


2.1 x 10- 2 


2.8 x 10~ 2 


8.0 x 10~ 3 


9.0 x 10- 1 


1.2 x 10" 2 


4.0 


4.8 x 10" 2 


4.0 x 10" 2 


1 x 10" 1 


2 x 10~ 2 


1.7 x 10" 2 


4.6 x 10~ 2 




|Re(^ 3 ) LR 




|Im(^ 3 ) LR 




0.25 


1.7 x 10" 2 


3.7 x 10" 2 


1.6 x 10" 2 


3.6 x 10" 2 


8.4 x 10" 2 


3.6 x 10~ 2 


1.0 


3.0 x 10~ 2 


3.3 x 10~ 2 


4.5 x 10~ 2 


6.6 x 10~ 2 


7.4 x 10~ 2 


1.0 x 10" 1 


4.0 


7.5 x 10" 2 


6.4 x 10" 2 


1.7 x 10" 1 


1.7 x 10" 1 


1.4 x 10- 1 


3.9 x 10- 1 




|Re(^ 3 ) LR 


= M^)rl 




|Im(^ 3 ) LR 


= |Im(^ 3 ) RL 




0.25 


1.4 x 10" 2 


5.9 x 10" 2 




2.3 x 10" 2 


4.4 x 10" 1 




1.0 


2.6 x 10~ 2 


5.2 x 10~ 2 




9.0 x 10~ 3 


2.3 x 10~ 2 




4.0 


6.5 x 10~ 2 


1.0 x 10" 1 




2.3 x 10~ 2 


4.4 x 10~ 2 





Table 1: Maximum allowed values for |Re(<5f 3 )AB| and |Im(5f 3 )AB|, with A, B = (L,R). A new 
parameter y is introduced, with y = m^/rn^. The definition of other parameters and their values 
follow Ref. [3]. 
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5.2 SUSY flavor change in Qq model 

In the Q 6 model the mass matrices of squarks in the flavor basis can be written as 

/ 1 \ 



10 , (5.98) 
\ y J 

q can be u or d, A can be L or R. Making the same unitary transformation on the squark fields 
as the ones on the quarks which diagonalize the quark mass matrices, we find the mass matrices of 
squarks in the SUSY basis (where the gluino coupling matrix is identity in the flavor space) to be: 







( 





°\ 1 


rh d )l L = T d P;(m d f LL P d O d = m\ L 


I+{yf-l)0 T d P* d 








1 P d O d 






\o 








m )lL 



( 5.43 • 10~ 5 1.31 • 10" 4 -0.0074 



(5.99) 



1.31 • lO" 4 3.17 • 10 4 -0.0178 
\ -0.0074 -0.0178 0.9996 

Note that this matrix is real, a consequence of the phase factorization of the fermion mass matrix. 
Similarly, 

/ 0.0367 -0.1339 0.1318 \ 



m l R 



i + (y d R -i) 



-0.1339 0.4891 -0.4816 
\ 0.1319 -0.4816 0.4742 / 



(5.100) 



{rhu) 2 LL = m\ L 



/3.85-10^ 6 7.74 -10" 5 -0.0020 



i + (vt-i) 



7.74 • 10~ 5 
-0.0020 



0.0016 
-0.0394 



-0.0394 
0.9984 



(5.101) 



2 

RR 



m 



uR 



( 1.46- 10^ 5 2.94- 10~ 4 0.0038 



i + (y u R -i) 



\ 



2.94- 10~ 4 
0.0038 



0.0059 -0.0768 
-0.0768 0.9941 



(5.102) 



K° — K° mixing via squark-gluino loops have several contributions. The most stringent limit 
arises from the (LL) — (RR) mixing, which requires [2] 

\(y d -i)\ 



< 0.23 



m 



500 GeV 



(5.103) 



(0.51 + 0.49 y d y/ 4 

Here we have assumed yf = y R = y d , and took the gluino mass to be equal to the first two 
family squrak mass. For first two family squark mass of 500 GeV, this translates to the limit 
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0.77 <y d < 1.24. For 1 TeV squarks, this limit is relaxed to 0.58 < y d < 1.48. We see that for y d 
order one, the most stringent limit on squark mediated FCNC is satisfied. 

The Qq model also generates significant (RR) (RR) contributions to the K° — K° mixing. We 
find 

0.68 <y d < 1.37 (5.104) 

for squark and gluino mass of 500 GeV. This constraint is also easily satisfied in the model. 
In the system, the analogous constraints are (from the (LL) (RR) operator) 



\(y° 



v. 



(0.53 + 0.47 y d )V* 



< 0.69 



m 



500 GeV 



(5.105) 



This limit leads to 0.48 < y d < 1.85 for squark-gluino mass of 500 GeV. The (RR)(RR) squark 
mixing gives no constraint from the Bd system. Similarly, there are no constraints arising from the 
D° system, nor from other type of operators in the model. 

In the leptonic sector, we find the (LL) slepton mixing (which is the same for the (RR) slepton 
mixing) to be 



[rhe) 



2 

LL 



m 



cL 



(vl 



I 2.93 ■ 10~ 4 
-4.02 ■ 10~ 3 
\ -0.0167 



-4.02 • 10~ 3 
0.0550 
0.2280 



-0.0167 \ 
0.2280 
0.9447 / 



(5.106) 



There are stringent constraints on the mixing parameter ((<5 6 )ll)i2 from the decay [i — > ej [23]. 
On the face of it, the mixing presented above would appear to be in mild conflict with data by a 
factor of few. However, since such a constraint is very week for the ((<5 e )i?_R)i2 mixing, we point 
out that the flexibility in the lepton sector mass matrix can be used to make the (LL) contribution 
small in exchange for larger (RR) contributions. That is, assume B -C B' in Eq. ( 11. 4ft . 



5.3 Left— Right squark mixing and a solution to the EDM problem 

So far we have ignored SUSY flavor violation arising from the left-right squark mixings. It turns out 
that these operators do not give significant contributions to meson-antimeson mixings, since such 
mixings have fermion chirality suppression. However, these mixings can generate new contributions 
to the neutron (and electron) electric dipole moments. Here we analyze constraints from the EDM 
and suggest a simple solution to the SUSY EDM problem. 

First, as shown in Ref. [8], the trilinear A-term induced phases align with the phases of the 
fermion mass matrices, even without assuming proportionality of the A-terms with the respective 
Yukawa couplings. This feature arises due to the phase factorization of the fermion mass matrix. 
Left-right squark mixings also receive contributions from the superpotential /z-terms. We derive 
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the mass matrix for the down squark sector to be: 



>d)l R = F{ 



d* 



( 5 d \ 

S d 
V a d J 



+ V2Ff 



/ \ 

o o p d 
\ o P' d o J 



(5.107) 



with 



F 1 = fJ, 3 v u3 + iii 3 v ul , 



F d 



-«ul- 



(5.108) 



After the unitary transformations to the left and the right squarks, corresponding to case (1), 
we have the (LR) mixing matrix in the flavor basis as 



rh d y LR = O d P d (m^ LR P d cO dc = F{ 



+V2Ffe i ^ 



-1.75 • 10" 4 4.14 • 10" 4 3.09 • 10~ 5 
-4.46 • 10" 4 3.84 • 10" 4 -5.45 • 10~ 4 
0.0078 -0.0286 0.0282 

4.53 • KT 5 -1.66 • 10 4 -1.24 • 10~ 5 \ 
1.79 • 10~ 4 -6.52 -10' 4 2.18 • 10~ 4 
-0.0031 0.0115 0.0125 / 



(5.109) 



with 



/ -1.62 • 10~ 5 2.18 -lO" 4 -0.0014 \ 



( m u) 2 LR — O^P u (mu) 2 LR P u cO u c — F{ 



+V2F?*e i A6u 



-2.18 -lO" 4 0.0022 -0.0283 
\ 0.0027 -0.0554 0.7168 / 

/ 3.24 • 10~ 5 -6.53 • 10- 4 0.0042 \ 

6.53 -10- 4 -0.0131 0.0848 . 
\ -0.0082 0.1661 0.0162 J 



(5.110) 



F± = n 3 v d3 + H3iv dl , 



F2 = VisVds + 



A*i + A*i2 



-Vdl- 



(5.111) 



(rheY LR = OiP e {m,)i R P e .O eC = F{ 



d* 



/ -1.10 • 10 5 2.85 • 10~ 4 6.72 • 10~ 4 

-2.75 • 10~ 4 0.0039 0.0257 
\ -4.89 • 10~ 5 0.0047 0.0313 

/2.93-10- 6 -1.14 • 10 4 -2.69 • 10" 4 \ 



+V2F" 



i A6» d 



1.10 • 10" 4 

y i.96 • 10- 5 



-0.0043 
-0.0019 



-0.0103 
0.0093 



(5.112) 
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Corresponding to case (2) these matrices are: 



\ m d)LR 



T d P d {m~)l R P d cO dc 



+V2Ffe* A8d 



( -2.25 • KT 4 
-5.73 • 10~ 4 
V 0.0100 

/ 5.20- 10~ 5 -1.90 



V 



2.05 • 10~ 4 
-0.0036 



5.32 • 10~ 4 
4.93 • 10~ 4 
-0.0368 

10~ 4 -1.42 



3.97 -10- 5 \ 
-7.00 • 10- 4 
0.0362 

• io- 5 \ 



-7.48 • 10~ 4 
0.0131 



2.50 • 10" 4 
0.0144 



7 



(5.113) 



( m u) 2 LR — O^P u (mu) 2 LR P u cO u c — F"* 



/ 



V 



-2.11 • 10~ 5 
-2.83 • 10~ 4 
0.0036 



2.83 • 10~ 4 
0.0028 
-0.0720 



-0.0018 \ 
-0.0368 
0.9319 



/1.41-10 -5 -2.83 -IO" 4 0.0018 \ 
-0.0057 0.368 
0.0720 0.0070 / 



2.83- 10" 4 
-0.0036 



(5.114) 



2 

LR 



mi 



) 2 LRPe c O e c — 



\ 



-1.41 • IO" 5 
-3.53- 10~ 4 
-6.28- IO" 5 



3.66 • IO" 4 8.62 • IO" 4 \ 
0.0050 0.0330 
0.0061 0.0402 



/ 



+V2Ffe i A8d 



/3.36-10- 6 -1.31 -10" 4 -3.08-10~ 4 \ 
1.26 • 10 4 

\2.24-10~ 5 -0.0022 0.0106 / 



-0.0049 
-0.0022 



-0.0118 
0.0106 



(5.115) 



Note that these matrices are in general complex, since F^' d are complex because of the spontaneously 
induced phases of the VEVs. This means that these matrices will contribute to neutron and electron 
EDM. Since these complex coefficients are proportional to fiv/m 2 , we find a simple solution to the 
SUSY EDM problem: Let the \x terms be of order 100 GeV, in which case one finds a suppression 
factor of 10~ 2 for the effective phase that enters the EDM expression. With this suppression factor, 
from the (1, 1) elements of these (LR) mixing matrices, we see that neutron and electron EDM 
constraints can be satisfied, even with the spontaneously induced phases in the VEVs being of 
order one. 

The proposed solution to the SUSY EDM problem has direct experimental consequences for 
LHC. We predict that the Higgsinos should be light, and three such pairs of doublet Higgsinos 
should be observable at the LHC. Their scalar partners, however, are inaccessible, since their 
masses lie in the few TeV range. 



30 



6 Conclusion 



In conclusion, we have presented a detailed analysis of the Higgs potential involving three pairs of 
Higgs doublets in a Q$ model of flavor. This class of models are motivated on two grounds: They 
lead to reduced number of parameters in the fermionic sector, and they can be helpful in alleviating 
the flavor changing problems of generic SUSY models. 

We have shown that tree-level Higgs boson induced FCNC are within experimental limits, even 
for the most stringent K° — K° mixing amplitude. The Higgs boson masses must lie in the TeV 
range. New sources of CP violation in meson mixing are highly suppressed. We have also shown 
the consistency of the model with SUSY flavor violation. A simple solution to the SUSY EDM 
problem is suggested, which requires light Higgsinos. 
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